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In this paper, we study phase diagram of a system of two-component hard-core bosons with 
nearest-neighbor (NN) pseudo-spin antiferromagnetic (AF) interactions in a stacked triangu- 
lar lattice. Hamiltonian of the system contains three parameters one of which is the hopping 
amplitude t between NN sites, and the other two are the NN pseudo-spin exchange interac- 
tion / and the one that measures anisotropy of pseudo-spin interactions. Physical states of 
the system are restricted to the ones without the doubly-occupied state at each site. We inves- 
tigate the system by means of the Monte-Carlo simulations and clarify the low-temperature 
phase diagram. In particular, we are interested in how the competing orders, i.e., AF order 
and superfluidity, are realized, and also whether supersolid forms as a result of hole doping 
into the state of the V3 x V3 pseudo-spin pattern with the 120° structure. 
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1. Introduction 

Study of systems in which competing orders coexist has been one of the most interesting 
topics in condensed matter physics. Among them, recent experiments on 4 He under pressure 
have led to renewed interest of supersolid. 1} Also the search for a lattice supersolid has been 
motivated by the realization of optical lattices in ultracold atomic systems. As dimensionality 
and interactions between particles are highly controllable and also there are no effects of im- 
purities, cold atomic system in an optical lattice is sometimes regarded as a "final simulator" 
for quantum many-body sy stems. 2) Numerical studies of hard-core bosons on a triangular 
lattice find that stable supersolid states form on doping of holes away from a 1 / 3 -filled solid 
(2/3-filled solid) with the V3 x V3 pattern. 3, 4) This provides example of old idea 5) that a 
finite density of defects in the solid (vacancies or interstitials) may Bose condense and form 
a superfluid in the existing crystalline background. In this paper we shall pursue the possibil- 
ity of realization of "vacancy condensation" phenomenon 6) in boson systems in a triangular 
lattice in which a frustration exists. 

The model that we study in this paper is a system of two-component hard-core bosons 
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with nearest-neighbor (NN) pseudo-spin antiferromagnetic (AF) interactions in a stacked tri- 
angular lattice. The particle number at each site is less than unity, and therefore the model is 
sometimes called bosonic t-J model, 7) i.e., a bosonic counterpart of the fermionic t-J model 
for the strongly correlated electron systems. The bosonic t-J model on a square lattice was 
originally studied as an effective model of the Bose-Hubbard model with strong on-site repul- 
sions, and its phase diagram was clarified. In this paper, we shall study the mode on a stacked 
triangular lattice. 

This paper is organized as follows. In Sec. 2, we shall introduce the bosonic t-J model and 
explain its path-integral formalism for the numerical investigation. In order to incorporate 
the local constraint faithfully, we introduce a slave-particle representation. In Sec. 3, we shall 
show results of the numerical study including phase diagrams, various correlation functions, 
and snapshots. We find that the system has a rich phase structure. Physical meanings of the re- 
sults is discussed by using, e.g., a mean-field theory like approximation. Section 4 is devoted 
for conclusion and discussion. 

2. Model and path-integral formalism 

The mode that we study in this paper is the bosonic t-J model 7) in a three-dimensional 
(3D) stacked triangular lattice whose Hamiltonian is given as follows, 

+/£(S?S? + S?s;), (1) 

(ij) 

where a] and b] are hard-core boson creation operators at site /, pseudo-spin operator 
Si = \b\&Bi with Bi = (a,-, b t )\ & are the Pauli spin matrices, and (/, j) denotes the nearest- 
neighbor (NN) sites in the 3D stacked triangular lattice. Physical Hilbert space of the system 
consists of states with total particle number at each site less than unity (the local constraint: 
a]cii + b]bi < 1). As we consider the 3D system, there exist finite-temperature (T) phase tran- 
sitions in addition to "quantum phase transition", which takes place as the parameters in H$ 
are varied. The system H$ might be derived as an effective model of Bose-Hubbard model 
that describes a cold atom system in an optical lattice. 8,9) In this case, / and J z are related to 
the intra-species and inter-species on-site repulsions, and / < 0. In this paper we study the 
bosonic t-J model (1) with general interests and treat / and J z as free parameters. We mostly 
consider the case J z , J > 0, i.e., the frustrated case because it is very interesting to see how 
the frustrated pseudo-spin state evolves as holes are doped. Relation between the geometrical 
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frustration and the formation of supersolid was discussed previously for the hard-core bosons 
in a triangular lattice. 3,4) 

In addition to its own theoretical interest, another motivation for studying the model 
(1) with the antiferromagnetic (AF) coupling is provided by its relation to the fermionic t- 
J model. Recently, the fermionic t-J model was studied by mapping it to a kind of bosonic t-J 
model by using a Chern-Simons gauge field. 10) The resultant bosonic model is closely related 
to the present model expressed by using the slave-particle expression given by Eq.(2) later on. 
In Ref., 10) the model was studied by a mean-field type approximation. Then numerical study 
of the model (1) and to obtain its "exact phase diagram" are important for understanding the 
phase diagram of the fermionic t-J model. Effects of the Chern-Simons gauge field will be 
discussed later on by comparing the obtained results in this paper and the phase diargam of 
the fermionic t-J mode on the triangular lattice. n) Furthermore very recently, it was shown 
by the numerical link-cluster expansion that the Fermi-Hubbard and Bose-Hubbard models 
on the square lattice exhibit very similar behavior in the strong-coupling region. 12) Therefore 
we expect that the results of the present study gives an important insight into phase diagram 
of the fermionic t-J model on the triangular lattice. 

In order to incorporate the local constraint faithfully, we use the following slave-particle 
representation, 

cii = 4>](pi U bi = <fi(p2i, (2) 

(<f>]fa + ip\xp Xi + ip\xp2i ~ l)lphys> = 0, (3) 

where fa is a boson operator that annihilates hole at site /, whereas (p^ (<x = 1,2) are bosons 
that represent the pseudo-spin degrees of freedom. Iphys) is the physical state of the slave- 
particle Hilbert space. Then the partition function Z at temperature T is given as follows in 
the path-integral methods, 

Z = J [D^D^D^e-^, (4) 

where /? = l/(k B T), H tJ is obtained by substituting the slave-particles representation (2) into 
Eq.(l), and the path integral is performed with satisfying the slave-particle constraint (3). The 
original path-integral representation of the partition function contains terms like aid T a i9 where 
at is the complex number corresponding to a] and r is the imaginary time. As we discussed 
in the previous papers 13) and also showed explicitly by the numerical studies on certain mod- 
els, 4, 1415) effect of non-zero Matsubara- frequency modes in the 3D system at finite tempera- 
ture is mostly the renormalization of the critical temperature and then the partition function 
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in Eq.(4) is a good approximation for studying phase diagram at finite-r. Furthermore, as the 
system in the 3D stacked lattice can be in a sense regarded as a sequence of the 2D system, 
its low-T phase diagram is closely related to that of the 2D system at T = 0. Therefore, it 
is expected that the low-T phase diagram of the 3D system in the stacked lattice obtained 
from Z in Eq.(4) is quite similar to that of the 2D system at T = as verified in the previ- 
ously studied cases. 4, 1415) In other words, the spatial third direction perpendicular to the 2D 
lattices plays a role similar to the imaginary-time direction. More detailed discussion on the 
model (1) at T = will be given in a future publication. Furthermore, it should be remarked 
here that a closely related discussion on the appearance of a Lorentz-invariant critical theory 
was given for the quantum phase model and the existence of Higgs mode was predicted in a 
two-dimensional superfluid on an optical lattice. 16) 

3. Numerical study 

We employ both the grand-canonical and canonical ensemble for the practical calculation. 
In the grand-canonical ensemble, the chemical potential term like /u (fit fa is added to H^. 
On the other hand in the canonical ensemble, the path integral in Eq.(4) is evaluated by 
means of the Monte-Carlo simulations with keeping the average density of holes fixed. To 
show results of numerical study, it is convenient to introduce the following dimensionless 
parameters, cj = fiJ, c Jz = fiJ z , c t = fit and a = J/J z . Therefore large cj (c Jz ) and/or c t 
corresponds to low-T region. 

To study the phase diagram, we calculate the internal energy E and the specific heat C 
defined as 

E=^(H tJ \ C = ^((H tJ - E)\ (5) 

where N = L 3 with the linear system size L. We performed calculation up to L = 30, 
and for the simulations, we employ the standard Monte-Carlo Metropolis algorithm with 
local update. 17) The typical sweeps for measurement is (30000 ~ 50000) x (10 samples), 
and the acceptance ratio is 40% ~ 50%. Errors are estimated from 10 samples with the 
jackknife methods. We also calculate the hole density p = (<f>]fa) 9 and correlation functions 
G xy (iJ) = <(S?S* + S y t S])) 9 G z (iJ) = (S'S)), and G s (iJ) = <£■ ■ Sj). It easily verified 
S 2 t = {a]ai + b]bi) 2 , and then the magnitude of the pseudo-spin is decreased by the doping 
of hole. On the other hand, boson correlation function G B (/, j) = {BjBj) is used to see if a 
superfluid forms. 
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Fig. 1. Phase diagram for anisotropic AF Heisenberg model in a stacked triangular lattice. Vertical axis plots 
cj z . At low T, there are two phases with a long-range order of 120° pattern, which are separated with each 
other by the line a = 1 (J = aJ z ). Phase transitions from paramagnetic to ordered states are of second order, 
whereas that of a = 1 is called morphotropic phase boundary. 19) System size L = 24 with the periodic boundary 
condition. 
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Fig. 2. E and C as a function of a for cj z = 1.5. At a = 1, C does not have a peak but exhibits some peculiar 
behavior though E indicates the existence of the phase transition. This is a typical behavior of the morphotropic 
phase boundary. 



3.1 Anisotropic AF Heisenberg model on the stacked triangular lattice 

At zero hole density, the system (1) reduces to the anisotropic AF Heisenberg model. 
The same model in the 2D triangular lattice has been studied rather intensively. 18) To identify 
the phase boundary of the present 3D model, E and C were calculated by the Monte-Carlo 
simulations. The spin correlation functions were also calculated in order to identify each 
phase. We show the obtained phase diagram in Fig.l. As T is lowered, second-order phase 
transition from the paramagnetic phase to the spin-ordered states takes place. Specific heat 
C exhibits a sharp peak at the transition points indicating the existence of the second-order 
phase transition. Location of the phase transition points is determined by the peaks of C. In 
the low-r region, there are two phases separated with each other by the line a = 1 .In the study 
of ferroelectric materials, the corresponding line is called the morphotropic phase boundary, 
and it plays an important role. 19) We show the internal energy E and the specific heat C as a 
function of the anisotropy a in Fig. 2, which exhibit a "typical" behavior of the morphotropic 
phase boundary . 9) It should be remarked that the system has the S U (2) pseudo-spin symmetry 
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at a = 1, whereas U(l) x Z 2 otherwise. This fact is the origin of the peculiar behavior of E 
and C at a = 1 . The phase structure at low T, i.e., o, o z — > large, is essentially the same with 
that of the 2D system at T = 0, as it is expected from the above general consideration. We 
have also studied the correlations of spins in the vertical direction of the stacked triangular 
lattice, and found the strong AF long-range correlations in the two phased with the 120° spin 
order and short-range one in the paramagnetic phase in Fig.l . 
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Fig. 3. Phase diagram for bosonic t-J model in a stacked triangular lattice at low T,cj = 6.0 (grand-canonical 
ensemble), p is hole density. Dots denote location of the phase transition observed by the numerical study. All 
phase transitions are of first order. SF stands for superfluid and this phase also has FM long-range order. The 
phase p ^ is essentially the pure spin system with the V3 x V3 symmetry and the phase with p ^ 1 is the 
empty state. The phase with p ^ 1/3 is explained in Fig. 6. These three phases are all insulating. 




3.2 t-Jxy model 

In this subsection, we focus on the xy AF Heisenberg model by setting / > 0, J z = in 
Eq.(l) and study effect of the hole doping. The case of / = 0, J z > will be investigated 
in the subsequent subsection. For the ferromagnetic case / < and |/| > \J Z \, it is readily 
expected that the system has a similar phase diagram to that on the cubic lattice 9) as there 
are no frustrations in the system. One of the motivations to study the doped AF magnets on 
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Fig. 5. Spin correlation functions Gs(/, j), G xy (i, j) and G z (/, j) for c f = 2.0 and fin = 12.0 (left), fifi = 14.0 
(right). 
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Fig. 6. (Left) Snapshot for c t = 4.5 and Pji - 14.0. Holes are localized and a kind of AF configuration of 
pseudo-spin is realized. We show pseudo-spins in the S x - S y plane for S z component is negligibly small. 
Therefore the length of arrows indicates magnitude of pseudo-spins. (Right) Caricature of typical configuration 
obtained by MC simulation. 



the triangular lattice comes from the possible spin liquid and superconductivity suggested by 
Anderson, 20) and therefore we focus on the AF cases. 

We first investigate the system in the grand-canonical ensemble. Obtained phase diagram 
in the c t - ftfi plane for cj = 6.0 is shown in Fig. 3. For the region c t < 5.5, there exist three 
phases and they are separated by sharp first-oder phase transitions. See calculations of E and 
the averaged hole density p = ((f>J(f)i) in Fig.4, the both of which exhibit sharp discontinuities 
at fifi 13.5 and 14.7. From this observation, we judge the existence of the first-order phase 
transitions there. The location of the phase transition points is determined by E. 

It is obvious that the phase for fifi < 13.6 is nothing but the pure spin system of very 
low hole density that has the long-range order of the three- sublattice V3 x V3 pattern. See 
the correlation function for fifi = 12.0 in Fig. 5. On the other hand in the intermediate region 
13.6 < fin < 14 A, stable state with p = - is realized. Correlation functions in Fig. 5 indicate 
that the state is nothing but the one shown by the snapshot in Fig. 6. We verified this conclusion 
by calculating the hole density. For \i > 14 .4 the density of hole is almost unity and the empty 
state forms there. 
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Fig. 7. Spin (left) and boson (right) correlation functions for c r = 12.0 and /fy/ = 1 1.0 



As c t is increased, all the above three phases make a phase transition to a new phase via 
first-order phase transitions. E exhibits discontinuities at the phase boundary of the phase 
p ~ 0.5 ~ 0.7 in Fig. 3, indicating first-order phase transitions. Hole density of this phase is 
p ^ 0.5 - 0.7, and the calculated spin correlation indicates the existence of the ferromagnetic 
(FM) long-range orders (S?) ^ and also non- vanishing superfluidity (Bj) ± 0. See Fig.7. 
This result means that the system is composed of, say, a bosons and the Bose condensation 
of a boson takes place there. If we impose the condition that the total number of a boson and 
that of b boson are equal, the phase separation to a-rich region and Z?-rich region is expected 
to occur. This problem is under study and result will be published in a near future. 

In the grand-canonical ensemble, the most stable state in the system appears for each value 
of the chemical potential. Near the first-order phase transition, it is rather difficult to control 
the particle density by varying the value of the chemical potential. Then it is quite interesting 
to study the system in the canonical ensemble by keeping the hole density constant. In the 
present system, we focus on how the state of hole density, say, 40% evolves as the hopping 
parameter t is increased. 
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Fig. 8. Internal energy E as a function of c t for p = 0.4 and cj - 10.0. At c t ^ 13, there exists a first- 
order phase transition. (Inset) N(E) for c t = 13.2 exhibits the double-peak shape indicating the first-order phase 
transition. 



In Fig. 8, we show the internal energy E for cj 
the number of states N(E), which is defined as Z 



= 10.0 as a function of c t . We also measures 
= J dEN(E) txp(-fiE). The result indicates 
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that there exist a first-order phase transition at c t = 13.2. In order to understand the physical 
meaning of the phase transition, it is quite useful to see snapshots of the two phases separated 
by the phase transition. See Fig. 9. From the snapshot for c t = 5, it can be seen that the phase 
of p = \ survives and there is a void of very low particle density as a result of an excess 
of holes compared to p = |. On the other hand, the snapshot for c t = 15 shows that the 
phase separation takes place, i.e., the region of pure-spin phase with the V3 x V3 pattern 
and the region of the superfluid coexist, but they are immiscible. In the superfluid region, the 
boson has a nonvanishing expectation value (B t ) ± 0. The observation obtained through the 
snapshots is verified by the correlation functions G B (/, j) shown in Fig. 10. 

The above result shows that the phase separation takes place and the supersolid does not 
form. The reason why the first-order phase transition takes place and the phase- separated 
state forms in the present model is understood as follows. Bose condensation {B t ) ± natu- 
rally induces a spin order. In the mean-field approximation, the wave function ^ B c of Bose- 
condensed state with a coplaner spin order in the S x -S y plane is given as 

^Bccc]^[^fl; + ^; + c]io>, 

i 

where c is a positive number. Then (S*)/(S y .) = cot(^/ - 0,-). On the other hand, (a t ) = ce lTli 
and (bi) = ce l6i . Therefore if the supersolid with the spin 120° long-range order forms, the 
phase of the superfluid cannot be uniform. As a result, the lowest hopping-energy state of the 
Bose condensate cannot be realized. More precisely, the expectation value of the Hamiltonian 
in the state ^ B c is evaluated as 

(Hu)bc tc 2 (cos Atj + cos Ad) + Jcos(Arj - AO) 

where A77 etc are the phase differences between Bose condenses on adjacent sites. In order 
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Fig. 9. Snapshots for p = 0.4, c/ = 10.0 and c t = 5.0 (left), c t = 15.0 (center, right (oblique angle)). Left and 
center snapshots show S x - S y component of pseudo-spin. Right one is snapshot from an oblique angle, and 
length of arrows indicates magnitude of pseudo-spins. 
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Fig. 10. Boson correlation functions for p = 0.4 and c t = 15.0 along the two lines shown in Fig. 9. 
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Fig. 11. Phase diagram for c t < 0. There are four phases as in the case of c t > and they are separated by 
first-order phase transition lines. Physical meaning of each phase is explained in the text. 



to generate the Bose condensation, the parameter t has to exceed some critical value. The 
hopping term with the coefficient t prefers A77, A6 ~ 0, i.e., the Bose condensation tends to 
accompany a ferromagnetic order. Then as c t is increased, a first-oder phase transition takes 
place and the system tends to phase separate into the superfluid region of intermediate particle 
density and the pure-spin region with 120° spin order and p ^ 0. 

From the above discussion, it is interesting to study the case c t = fit < 0, which is 
sometimes called frustrated NN hopping. From the above consideration, one can expect that a 
state with both a non-colinear spin order and the superfluidity with a nonvanishing momentum 
(i.e., 7T > IA77I, |A#| > |) forms at sufficiently low T. Here it should be mentioned that the case 
c t = fit < can be realized in a rotating Bose gas system as rotation of optical lattice generates 
an effective magnetic field for bosons. 21) In the case in which magnetic flux penetrating each 
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Fig. 12. Spin (left) and boson (right) correlation functions for c t = -19.0 and /3/i = 1 1.0. The spin correlation 
function G xy (/, j) exhibits the V3 x V3 pattern. The boson correlation Gb(/, j) also shows a similar behavior. 
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triangular plaquette is exactly 7r, the frustrated NN hopping is realized. Furthermore in the 
fermionic t-J model, the model with t < describes the electron doped materials. As we 
mentioned before, the fermionic t-J model is related to its bosonic counterpart through a 
Chern-Simons gauge theory. Possible relation of the bosonic t-J model to the fermionic t-J 
model will be discussed later on. 

We numerically studied the t-J model with a negative t as in the previous case with t > 
0. We first show the obtained phase diagram for c t < in Fig.ll. There are four phases 
and they are separated by first-order phase transitions as in the previous cases. The phase 
transition lines are determined by the measurement of E. The phases with p ^ 0, |, 1 are 
essentially the same with the ones shown in Fig. 3. The new phase that appears for c t < -7 
is the expected to have both the co-planer long-range spin order and the superfluidity. To see 
it, we show the spin and boson correlation functions in Fig. 12. It is obvious that the pseudo- 
spin has the long-range order with V3 x V3 pattern and also the Bose condensation with a 
nonvanishing momentum forms there. Appearance of this phase comes from the fact that the 
Bose condensates {a t ) and {b t ) have different phases depending on the A, B and C sublattices, 
and these position-dependent condensations are enhanced by both the 120° structure of the 
spin oder and the negative hopping t < 0. This the reason why this phase is different from 
that for large positive t in Fig. 3. 

The correlation function G z (/, j) in Fig. 12 exhibits vanishing value and therefore the den- 
sity of atoms are uniform even for c t < -7. This means that supersolid does not form in the 
present system. In the following subsection, we shall study the t-J z model and show existence 
of the supersolid phase as a result of the competition of t and J z terms. As the t-J z model with 
J z < is directly derived from the Bose-Hubbard model with repulsions, the obtained phase 
diagram is expected to be verified by experiments of the two-component cold atom systems. 

3.3 t-J z model 

In the previous subsection, we studied the t-J^ model and clarified its phase diagram. 
There we found that there the supersolid does not form though the phase with the both the 
spin and boson long-range orders exists in some parameter region. In this subsection, we shall 
continue the numerical study on the bosonic t-J model and consider the case with J z > 0, / = 
0. The case J z > is realized in the two-component cold atom system in which the intra- 
species repulsion is larger than the inter-species one. 9) As S? = \{a]a t - b]bi), the existence 
of the long-range orders (S?) ^ 0, {a t ) ± {{b t ) •$> 0) induces a genuine supersolid. 

The obtained phase diagram by the MC simulations is shown in Fig. 13. Phase transition 
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Fig. 13. Phase diagram of t-J z model in the grand-canonical ensemble. There are five phases and they are 
separated by first-order phase transition lines. SS stands for supersolid. cj z = 10.0. 
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Fig. 14. Behavior of internal energy and hole density at phase boundaries of supersolid indicating phase 
transitions at c t ^ 12.0 and c t ^ 20.0, where J3/J, = 5. Arrows indicate the location of phase transitions. 



from the pure-spin state to the supersolid is of second order and the others are of first order. 
The locations of the phase boundary are determined by the measurement of E and C. As 
we expected, there is a parameter region in which the supersolid forms. The behavior of the 
internal energy and hole density are shown in Fig. 14 at the phase boundary of the supersolid. 
Correlation functions to be used for identify the supersolid are explicitly shown in Fig. 15. 
The correlation of the z-component of spin indicates there exists solid order, and the particle 
correlation shows the existence of a small but finite long-range order, i.e., superfluidity. 
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Fig. 15. Spin and particle correlation functions in supersolid. c t = 17 and fifi - 5. 
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4. Conclusion and discussion 

In this paper we studied phase diagram of the bosonic t-J model in the stacked triangular 
lattice. The model has a rich phase structure and we expect that some of them are observed by 
experiment on systems of two-component cold atomic gas and strongly correlated electron 
systems. 
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Fig. 16. Phase diagram of the hard-core boson system in the stacked triangular lattice. 4) In this figure, p 
denotes the particle density. V and t are parameters of the repulsion and hopping, and we set fiV = 50 (low-T 
region). There exists a supersolid state between two solid states. 
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Fig. 17. Density profiles of a, Z?-atoms and hole in the supersolid. cj z = 10.0, c t = 17.0 and/fy/ = 5.0. 



It is interesting to compare the obtained phase diagrams with that of the hard-core bo- 
son system in Ref., 4) which is shown in Fig. 16. In the model Hamiltonian, V represents the 
nearest-neighbor repulsive interaction and t is the hopping parameter. There are two solid 
states with the particle density ~ 1/3 and 2/3, and a supersolid state between them. The 
phase diagram in Fig. 16 respects the particle-hole symmetry of the model, whereas in the 
t-J model such a symmetry does not exist. In the t-J^ model studied in the paper, the super- 
solid does not form in similar parameter region and the phase separation takes place there 
instead, whereas in the t-J z model the supersolid appears as in the hard-core boson system. 
For t-J^ model with t > 0, the reason why the phase-separated state dominates was explained 
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in Sec. 3. 2. The state shown in Fig. 6 corresponds to the state with the particle density = 2/3 in 
the hard-core boson system. In the t-J model, the state with particle density =1/3 generally 
does not have spin long-range order nor solid (density-wave) order because of the low particle 
density and the non-existence of the particle-hole symmetry. Structure of the supersolid in the 
t-J z model is slightly different from that observed in the hard-core boson system. The particle 
density there is larger than 2/3, and one kind of atom, e.g., £>-atoms occupy 5-sublattice and 
a-atoms and holes form a superposed state in A and C-sublattices. See Fig. 17. Therefore for 
appearance of the supersolid in the t-J z model, spin degrees of freedom plays an essential 
role. On the other hand in the hard-core bosons on the triangular lattice, the supersolid can be 
interpreted as a superfluid of excess bosons (or holes) via a relay-type movement. 

It is interesting and useful to discuss the relation between the bosonic t-J model studied 
in this paper and the fermionic t-J model on the triangular lattice whose phase diagram was 
studied by high-temperature expansion in Ref. n) For the Hubbard model on the square lattice, 
very recent study by the numerical link-cluster expansion shows that the Fermi-Hubbard and 
Bose-Hubbard have a very similar physical properties for large on-site repulsions. 12) On the 
other hand in Ref., n) the authors studied the fermionic t-J model on the triangular lattice and 
concluded that the RVB state appears by hole doping into the state with the three- sublattice 
spin symmetry for t > and p = 0.2 - 0.6. By the high-temperature expansion, they found 
that the entropy decreases considerably at low temperature and peak of the spin susceptibility 
moves to the high-temperature region. These results indicate that the hole doping releases 
the frustration and enhances AF correlation though AF long-range order is not observed. 
In their discussion on the realized state, it was assumed that the low-temperature state is 
homogeneous and metallic. If this assumption is correct, the conclusion that the RVB forms at 
low-temperature seems plausible and acceptable. However the results obtained in this paper 
clearly offers another interpretation of their calculations, i.e., the state shown in Fig. 6 has 
obviously has lower entropy and stable AF correlation compared to the 120° spin state at 
p = 0. Holes are localized there and the state is inhomogeneous contrary to the assumption in 
Ref. n) We expect that a similar state to that in Fig. 6 is realized in the fermionic model, though 
the superfluid in Fig. 3 corresponds to a metallic phase of mobile holes in the fermionic t-J 
model. On the other hand for t < 0, it was found in Ref. n) that the entropy remains large until 
very low temperature and therefore the hole doping does not release the spin frustration at 
p = 0. This result is also consistent with the results shown, e.g., in Fig. 12. 

In the above discussion on the bosonic and fermionic t-J models on the triangular lattice, 
we directly compared the obtained results for the models, and concluded that there is a close 
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resemblance between two models. It is important and interesting to show the relationship 
between two model by an analytical method. To this end, discussion using the Chern-Simons 
gauge theory in Ref. 10) is useful. In the slave-particle representation (2), the Hamiltonian (1) 
is expressed as 

H t j = -t J] (<pli<f>i<f>)<P<rj + h.c.) 

(i,j),cr 

+J ^ (^(pl^lje^^cp^j^i + h.c.)PiP j9 (6) 

(i,j),(r,cr' 

Pi = l-(f>](f>i, (7) 

where & = 2(1) for <x = 1(2), e n = -621 = 1, £11 = ^22 = 0, and we have set J z = J 
for simplicity. On the other hand, the Hamiltonian of the fermionic t-J model in two spatial 
dimensions is given as follows by using the Chern-Simons gauge fields A*, and A*, defined on 
the link 0',;'), 10>22) 

(ij),cr 

+J J] (^^^/^^^^/ + h.c.)P i i > y, (8) 

(i,j),cr,cr' 

where crA = A (-A) for <x = 1 (2), and 

^ Ay = tt ^ ((f\(fi - ^ 2 ), mod 2tt, 

(i,j)eC inside C 

£ A* = 7T £ (0V - P)> mod 2tt, (9) 

(iJ)eC inside C 

where C is a closed loop and Ay = -A j7 is understood in the above summations in (9). 
In the homogeneous state, we can replace the RHS's of Eq.(9) by their mean values for 
investigating the ground-state phase diagram. In the case of the coplanar spin configuration 
as we considered in this paper, effects of At is small because (((p\<pi - (p\yi)) = 0. Similarly 
for a homogeneous hole distribution, the Chern-Simons gauge field Ay can be set Ay = 
in the mean-field level. Therefore for a homogeneous state, phase diagrams of H t j and 
are closely related with each other, though dispersion relation of excitations in the two 
systems are different by the local constraints Eq.(9), i.e., in the system H^, hopping of bosonic 
spinon and holon accompanies the Chern-Simons flux. 23) Anyway, more detailed study on the 
relation between the two models by using the Chern-Simons gauge theory is interesting and 
useful. 24) 
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